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ABSTRACT 

We analyse the role of a fluctuating a-effect in a£l dynamo models, and show that 
there is a mechanism for magnetic field generation, valid at large scale separation, 
deriving from the interaction of mean shear and a fluctuating a-effect. It is shown 
that this effect can act as a dynamo even in the absence of a mean a-effect, and that 
the timescale for dynamo waves is strongly affected by the presence of fluctuations. 
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1 INTRODUCTION 

The mean-field ansatz has for many years been an invalu- 
able tool in the construction of tractable dynamo models for 
the Sun, stars and planets. Extended treatments of the the- 



ory may be found in the monographs of Moffatt ( 1978 1 and 
Krause fc Radler] ( |1981[ ). In its traditional form, the theory 
distinguishes sharply between helical and non-helical turbu- 
lence. The latter case is presumed to imply that the statistics 
of the small-scale velocity field that ultimately produces the 
mean-field effects are symmetric under reflection; this leads 
to the vanishing of the a-effect term in the mean induction 
equation. Any mean inductive effects of a more general sort 
(for example, the "f2 x J" effect of Radler) are associated 
(like the diffusion term) with two spatial derivatives of the 
magnetic field, and are not, unlike the a-effect, guaranteed 
to lead to mean field growth on sufficiently large scales. 

Investigation of the dynamo properties of velocity fields 
which are non-helical on the average, however, has usually 
ignored the possible effect of fluctuation. It is well known 
that even non-rotating turbulence has large fluctuations of 
helicity, though the mean is essentially zero. Even if the sign 
of the helicity is well defined, there is no guarantee that the 
emf's generated from a mean field by the helical flow will 
not have large fluctuations. The paper of Brandenb urg fc| 
Sokoloff ( 2002 I , which made a systematic attempt to derive 
transport coefficients from a simulation of shear flow tur- 
bulence, noted the presence of significant fluctuations. This 
effect is even more graphically illustrated by the results of 
a recent numerical experiment on convection in a rotating 
layer ( Cattanco & Hughes 2006 ) . A weak uniform horizontal 
field is imposed on a fully developed convective flow in the 
layer, and the emf measured. It is found that although the 
helicity, when averaged over the upper part of the convec- 
tive layer, has a definite sign, the emf fluctuates wildly in 
magnitude and direction and if there is a non-zero average it 
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Figure 1. Time series showing the behaviour of the mean emf 
(averaged in the horizontal and over the top half of the layer) for 
the rotating convection simulation of Cattaneo &; Hughes 2006). 
The three plots show the measured emf in each of the three coordi- 
nate directions (z is vertical). The heavier line shows the running 
time-average. The imposed magnetic field is in the x direction. 



is very small compared with the rms fluctuations. Figure [T] 
from their paper, shows this clearly. 



It is of interest, therefore, to understand whether large 
fluctuations in the induced mean emf, on a timescale slower 
than the original averaging process but faster than the 
timescale for mean field evolution, can lead to significant 
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dynamo action. A number of authors have considered this 
problem from various points of view. An early investigation 



was made s ome time ago by Kraichnan ( 1976 1 (see also Mof- 
fatt (19781, Ch. 7), who examined an a dynamo, with no 
mean flow, in which the a-effect exhibited fluctuations with 
zero mean. The analysis, based on first order smoothing, led 
to a possible negative diffusion effect, giving rise to the possi- 
bility of mean field growth. In the case of negative diffusion, 
however, the consequent rapid growth of small scales would 
violate the scale separation on which the averaging process 
was based. Thus although very interesting, the results could 
not be said to be entirely conclusive. 

The Cattaneo and Hughes calculation does not pos- 
sess large scale coherent shear, and the situation is there- 
fore somewhat like the scenario of Kraichnan, which started 
from the supposition of an a 2 -dynamo model. Indeed Catta- 
neo and Hughes found no large-scale dynamo mode in their 
calculation, in spite of the large rms a. The situation ap- 
pears rather different when there is a coherent zonal shear, 
such as occurs in the Sun. In that case there is a powerful 
mechanism for creating zonal from meridional field, and this 
might make the effects of fluctuating helicity more impor- 
tant. 

There have been a number of attempts to address this 
situation. Vishniac & Brandenburg ( 1997 1 investigated a re- 



duced mean field model appropriate to galactic dynamos. 
This takes the form of an ODE and the fluctuating a ef- 
fect plays the role of multiplicative noise. They show nu- 
merically that growth can occur for large enough fluctu- 
ations. They also report on a numerical experiment on a 
spatially extended model, but do not attempt a systematic 
survey. These ideas were further investigated by Fedotov et 



al. (20001 from the point of view of stochastic calculus. An 



approach close to that of the present paper was adopted by 
Silant'ev (20001. After very lengthy and detailed analysis it 
was concluded that the effect of fluctuations was a term of 
a-effect type, derived from spatial variations of the fluctua- 
tions. 

The purpose of this paper is to revisit the problem us- 
ing a simpler formulation than that of Silant'ev . We show 
that the situation is not quite as he envisaged. There is in 
fact a mechanism that leads to dynamo action, even when 
the mean value of a vanishes, and the fluctuations have no 
spatial dependence. The mechanism is guaranteed to work 
for sufficiently large-scale mean fields, in contrast to Kraich- 
nan's mechanism. In view of the numerical demonstrations 
described above, which show that relying on a well defined 
a may be misleading, it would seem that this new mecha- 
nism may well be of great importance in understanding the 
dynamo properties of turbulent flows. 

We demonstrate the mechanism by starting with the 
standard aQ dynamo equations, and allow a to vary rapidly 
in time. A simple asymptotic theory yields the new ansatz. 
The results are then applied to a simplified one-dimensional 
dynamo model to show how dynamo action is enhanced by 
the new mechanism. Finally we give results from a nonlinear 
one dimensional dynamo simulation incorporating the new 
mechanism, and demonstrate the important effect that the 
fluctuations have on the frequency of activity cycles. The 
paper concludes with suggestions for future development. 



2 DERIVATION OF THE MEAN FIELD 
EQUATIONS 

In order to explain the mechanism we reduce the model 
to its simplest form. Consider an axisymmetric mean field 
B = B(r,8)e^ + B p [= V x {A(r,8)e <t> )] in spherical polar 
coordinates (r,6,cj>). The only mean flow is that of zonal 
shear with differential rotation f2(r, 6). Due to small scale 
fluctuations there is a zonal emf proportional to B, of the 
form a(r,8,t)B . Here and throughout the fluid is taken to 
have uniform magnetic diffusivity r\. The induction equation 
then takes the standard form for an aQ, dynamo, 

rl A 

aB + r]D 2 A, 



rsintfBp • Vft + i]D B, 



(1) 
(2) 



dt 
dB 
dt 

where D 2 = V 2 — 1/r 2 sin 2 8. It is now envisaged that a 
varies in time on a timescale r that is rapid compared with 
the evolution time of the mean field (this can always be 
achieved provided that the scale of the field is large enough). 
We then define a small parameter s and write 

d d _i d 



Qo + e ai(r), with (qi) = 0, (3) 



d_ _ 
dt dt +S cV 
and (■} denotes an average over the fast timescale. We also 
define B = B + eBi(r), A = A + Ai(t), where the de- 
pendence on r, 6, t has been suppressed. The leading order 
fluctuating parts of Equations (Tj2l then give 

{ ^-=a 1 B , ^r 1 =rsin#B lp - Vfi. (4) 
or or 

When this is solved for Ai , B\ we can calculate the order 1 
quantity (aiBi) that will appear in the average equation for 
Ao. This is conveniently given in terms of the function y(r), 
defined by dy/dr — ax, {"/) ~ 0; we find Bi p = V x yB^e^ 
and 

(aiBi) = -<75i T ) = -(-yrsineV x yBoe^} ■ X7Q. 
= -<7 2 }V(rsin6iBo) x ■ Vft 

1 o 

- o V(7 ) x rsin6>B e^ ■ Vfi 



(5) 



The mean field equation for A can now be obtained from 
Equation |TJ by averaging, to yield 

8A 
dt 



= (aiBi) + a a B + rjD 2 A , 



(6) 



where the averaged term is given by Equation (|5| . The mean 
version of Equation |2| is unaltered except for subscripts 
zero on B p , B. 

We can explain this effect physically as follows: Con- 
sider a zonal magnetic field varying linearly in latitude. Then 
the action of a spatially uniform but time dependent a-effect 
is to produce a spatially uniform radial field that also os- 
cillates. This is turn is acted upon by the zonal shear to 
produce an oscillating zonal field. This is acted on by the 
fluctuating a-effect and the resulting zonal emf has a non- 
zero time average. 

It is not difficult to incorporate the effects of spatial 
variation of the fluctuations into teh analysis. The results 
are not qualitatively affected by this generalization, provided 
that the effects of spatial diffusion are smaller on average 
than that of the time dependence. 

It will be shown in the next section that the new terms 
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in Equation |6| are capable of producing growing fields even 
when ao = 0. We can see immediately that both the new 
mean term in Equation Q and the flux stretching term 
in Equation |2]| have only one space derivative in Bo, Ao 
respectively. This means that at sufficiently long length scale 
these terms will dominate over the diffusion terms. 



3 A ONE-DIMENSIONAL DYNAMO WAVE 

Rather than conduct a fully two-dimensional numerical ex- 
periment, we can show how the new term leads to field 
growth by considering a one dimensional wave model of the 



type originally proposed by Parker (see, for example, Proc 



aB + V (A XX - I 2 A), 



(7) 



tor & Spiegel (199f )). Specifically, we consider A and B to 
be functions of x, t alone, and consider the model 

OA 

~ = Q'A x +ri{B xx -e 2 B), (8) 

Here a and SI' are independent of x, though a is a function 
of time in the manner discussed above. Applying exactly 
analogous analysis to this equation, we obtain the following 
simplified version of Equation Q (the subscripts zero on 
A, B have now been dropped) 



-G 2 fl'B x + a B + rj(B xx - i 2 B), 



(9) 



9) 



dA _ 
dt ~ 

where G 2 = {j 2 }- This equation can now be solved together 
with Equation jjjj, to find marginal solutions of the form 
A, B oc expi(kx + cut). Substituting into equations 
yields the dispersion relation 

(iu> + r/(k 2 + I 2 )) 2 = G 2 n' 2 k 2 + iQ'aok. (10) 

The imaginary part lu — fl' aok/2i](k 2 + £ 2 ) (so that the 
dynamo waves travel if Sl'ao ^ 0), while the real part can 
be written in dimensionless form. If we write k = £m, and 
define 



V 



fl'cto 



2 



G 2 Vl' 2 
rj 2 P 



(11) 



then w = r)l 2 T> 
by 



m/2(l + m 2 ) and T>, Q and m are related 



^->2 2 

Q m 



v- 



4(1 + ', 



(1 + m 2 



(12) 



When Q — we have the usual aQ, dynamo. Dynamo ac- 
tion is possible when \T>\ 32/3-\/3, with equality when 
m — l/v3. Even when T> = 0, we can find dynamo action 
when \ Q\ ^ l/v2, with equality when m = 1. Although in 
this case lo — so that the wave does not travel, the ex- 
ample shows that the fluctuations in a alone can lead to 
growing magnetic fields, even when there is no mean emf. 
The stability boundary is shown in Figure [2] 



4 A SIMPLE MODEL OF THE SOLAR 
ACTIVITY CYCLE 

For dynamical reasons it has long been accepted that any 
a-effect should be antisymmetric about the solar equator. 
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Figure 2. Relation between T> 2 and Q 2 for marginal stability 
according to Equation \12\ . Each straight line is for a constant 
value of m, between m = 0.36 and m = 1.2. Dynamo action is 
possible above the envelope of the lines. 



The new term is related to the mean squared value of a 
(actually to the mean square of its time integral) , and thus 
on average (the statistics of a being odd, and so those of a 2 
even about the equator), one should take it as symmetric 
about the equator. This then, as one might hope and ex- 
pect, preserves all the parities of conventional models. To 
investigate cyclical behaviour we solve numerically a model 
consisting of Equation (JsJ) with f2' = rj = I = 1 together 
with the following version of Equation Q: 

rB x + dsin(2nx/l)B 



dA 



l + B 2 



+ B XX -B), (13) 



where d, r are positive constants. We have adopted the 
simplest possible quenching formula so as to ensure equi- 
libration of the mean fields. These equations are solved in 
< x < I, where A, B are made to vanish at x = 0, 1. It is 
easy to find time-periodic solutions of dipole parity with B 
antisymmetric about x = 1/2 . In Figure [3] are shown space- 
time plots of B(x,t) for / = 10, d = 30 and various values 
of r. It can be seen that as r increases and fluctuations be- 
come more important the period of the cycle rises rapidly, 
and eventually the solution becomes steady. 



5 DISCUSSION 

In this short paper it is shown that the interaction of shear 
with fluctuating mean field effects leads to a dynamo mech- 
anism that may well prove more powerful than the usual 
a-effect in astrophysically interesting situations, though the 
cyclical behaviour of any model will continue to depend on 
the usual a-Q interaction. While the derivation of the mean 
field equations has necessarily involved particular scalings 
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Figure 3. Space time plots of B(x,t) for solutions of Equations 
{§} ,{13} for I = 10, d = 30 and four different values of r. For 
r = we have the usual a — f2 dynamo. The cycle time increases 
with r (note the different timescale in the third image) for r large 
enough the solution becomes steady as shown in the final picture. 



for the amplitudes of the fluctuations, it should be empha- 
sized that these are in some sense arbitrary as the mech- 
anism can be guaranteed to give dynamo action at long 
enough scales, however small the effect. It is now necessary 
to integrate more realistic models in spherical geometry, and 
to derive the results for the mean field directly from the ve- 
locity fields; these, and the effects of spatial fluctuation, are 
the subject of work in progress. 
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